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Abstract 

Large values for the mass-to-light ratio (T) in self- gravitating 
systems is one of the most important evidences of dark matter. We 
propose a expression for the mass-to-light ratio in spherical sys- 
tems using MOND. Results for the COMA cluster reveal that a mod- 
ification of the gravity, as proposed by MOND, can reduce signifi- 
cantly this value. 
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1 Introduction 

The standard model for cosmology states that 95% of all content of the 
universe belongs to a dark sector, not yet detected, composed by a 
unknown form of energy called dark energy and some kind of non- 
baryonlc matter or dark matter, while baryogenesis implies that only 
5% of the universe is made up of atoms. In general terms, dark mat- 
ter is the unknown fraction that appears in agglomerated matter, like 
galaxies and clusters of galaxies, having a effective null pressure and 
suffering consequently the clustering process. Dark energy is a smooth 
component of the matter content of the Universe which must have neg- 
ative pressure, inducing the accelerated expansion of the Universe. A 
review of the observational evidences and theoretical explanations for 
dark matter and dark energy can be found in references (TJ |2j [3). 

The first evidence on the existence of dark matter came from the 
analysis of the dynamics of the Coma cluster of galaxies, in the thirties 
l4l . Using the virial theorem and the hypothesis that the Coma cluster 
is in dynamical equilibrium, it was found that there should have much 
more matter than it could be deduced using the usual mass-to-light 
ratio for galaxies. Later, the study of the rotation curve of spiral galaxies 
has shown that even at the level of galaxies the major part of matter 
does not emit any kind of radiation. In galaxies, the mass-to-light ratio 
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using the rotation curves can be as high as 30, in solar unities, while 
this number is of the order of 180 for the Coma cluster (5). 

In general, the existence of dark matter is proposed in order to ex- 
plain anomalies observed in the dynamics of self- gravitating systems, 
like galaxies and clusters of galaxies. An alternative to dark matter is 
to admit that the fundamental laws of physics are not the usual one 
from a certain scale on. This can be achieved, for example, by chang- 
ing the gravitational potential. Another possibility, at newtonian level, 
is to modify the Newton's dynamic laws. One interesting proposal in 
this sense is the MOND theory |6|, which is based on a changing of the 
Newton's second law, such that it takes the form 



where ^i(x) is a function such that ^i(x) sa 1 for x >> 1 and fi(x) sa x for 
x « 1. Hence, a defines a critical acceleration below which the usual 
Newton's second law is not valid anymore. This phenomenological pro- 
posal has been very succesfull in explaining the rotation curve of the 
spiral galaxies, a taking a universal value around a ~ 10 -8 cm. s~ 2 0. 
Some interpolation formula between the "newtonian" regime and the 
"MOND" regime for n(x) have been proposed in the literature (8). 

One of the main criticism with respect to the MOND theory is that 
it is non-relativistic, and it can not be deduced from a Lagrangian. 
More recently, a relativistic version of the theory has been proposed, 
involving, beside the usual geometric and matter sector, a vector field 
and a scalar field 1 9 1 . But, it remains some controversy if MOND can 
succesfully explain the lensing effect in clusters of galaxies 1 10 1. 

The proposal of the present work is to go back to the dynamics of 
cluster of galaxies, in the realm of the MOND theory. There has been 
some studies in this sense using X-ray data 1111 . The purpose of the 
present work is more modest: to generalise the virial theorem using 
a suitable interpolation expression for fi(x) and from it to deduce the 
mass-to-light ratio. To be specific, we work with a given set of obser- 
vational data concerning the Coma cluster. It will be shown that, even 
if a closed expression for the generalised virial formula is not possible, 
we can obtain the mass-to-light ratio that is compatible with ordinary 
estimations without dark matter: values as low 10 can be easily ob- 
tained. The paper is organised as it follows. In next section, we present 
the general expressions for virial theorem. In section 3 the observa- 
tional data for the Coma cluster are discussed and a estimation for 
the mass-to-light ratio is obtained. In the last section we present our 
conclusions. 
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Figure 1: Behaviour of the function fi(x). The newtonian limit is 
reached when fx -> 1. 



2 Method 

The general form of the MOND theory is given by 

m(J,(x)a = —G —^- f , (2) 

where (j, is a function, to be specified, which has the following proper- 
ties: fj,(x) ~ 1, for x » 1, and fi(x) ~ a; for x << 1. Since x — a/a , 
where ao is a constant with dimension of acceleration, the MOND the- 
ory implies the introduction of a new fundamental constant: a critical 
acceleration ao which specifies when the newtonian mechanics is valid. 
Expression (2) describes the motion of a point mass m under the in- 
fluence of the gravitational field of a spherical large mass distribution 
with total mass M. The function fi(x) is arbitrary, excepted by the fact 
it must obey the assymptotic regimes described before. We will adopt 
the following form for this function: 



(3) 



with 



x = — . (4) 

a 

For technical reasons, it is convenient to rewritte equation (2) as 

X= -G -r-5 — f. (5) 

H[x)r z a Q 

Some simple manipulations lead to the expression 



GM r 2 a „ 

aM = -—y 1+ GM r ■ (6) 
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With the identification y/GM/ao = r c , where r c is a "critical distance" be- 
yond which the MOND regime becomes effective, equation (6) becomes 




am = 5-Wl+-5-r . (7) 



It must be stressed that the MOND critical parameter is the accelera- 
tion scale ao. However, with the formal definition of a distance scale, 
the MOND theory is recast as a (formal) modification of the gravitational 
potential. The new form of the gravitational potential is dictated by the 
choice of the function (j,(x). We will use this approach from here on, but 
the results does not depend on it: it is just a question of description of 
the dynamics. Using such approach the potential corresponding to the 
function fi(x) chosen is 



$m (r) = -GM< 



1 + sMT 1 ^ 



y (8) 



The potential has the usual newtonian form in the limit r « r c , while it 
becomes logarithmic when r >> r c , what corresponds to the full MOND 
regime. 

The first indirect detection of dark matter came from the application 
of the virial theorem, which is valid in the Newtonian dynamics, to the 
Coma cluster |4|. This theorem establishes that in a self-gravitating 
system, in dynamical equilibrium, the potential energy is minus two 
times the kinetic energy. Measuring the kinetic energy, what can be 
done through a spectroscopic analysis of the light coming from a dis- 
tant system like a cluster of galaxies, the potential energy, and hence 
the mass of the system, can be deduced. But, such simple relation 
between the potential and kinetic energy can only be possible in the 
case of a newtonian theory. For a different gravitational potential (or a 
different dynamic law), the virial theorem in its usual form is not valid. 
The general relation for the virial theorem assumes the form 

K lyyfltMV 



dr 

where it was computed the contribution from all i-components of the 
system with mass Introducing (8) in the above relation the kinetic 
energy of a self-graviting system orbiting in the full MOND regime can 
be obtained as 

/<=«R, (10) 
2r y r A c 

where m is the mass of a tiny shell at radius r, and M (r) is the mass 
enclosed by this shell. 
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If we know the kinetic energy of a spherical system it is possible 
to evaluate the mass-to-light ratio T under the assumption that T is 
independent of radius. The kinetic energy can be associated with the 
surface luminosity I and the line-of-sight velocity dispersion a by 

K = TJ, (11) 

where J is a integral defined as 1 12 1 

J = 3tt / I(R)a 2 (R)RdR. (12) 
Jo 

Moreover, a symmetric mass distribution should be described by some 
kind of density profile. The mean density may be associated with the 
observational data and can be written as 1121 

T f°° dl dR 

The above definition enable us to calculate the expression for the mass- 
to-light ratio with MOND. Inserting fl3] twice in the MOND expression 
for the kinetic energy 1101 and combining this result with 11 U it is 
possible to determine the mass-to-light ratio using only the velocity 
dispersion and the surface luminosity data: 

r M = ~. (14) 

j 

In above relation J is the integral define before and J is defined as 

J=-16G [ t{r)p{r)rJl+ r -^dr, (15) 



jo y c 

where t(r) and p(r) are also auxiliar integral expressions that depend 
on / and a 

poo it jn pT 

w-lmvvT? *4 wV ' 1161 

In order to calculate T for a self-gravitating system we need the surface 
luminosity profile I(R) and the dispersion velocity profile a(R). Note 
that the velocity dispersion appears only in term J. Hence, J represents 
the kinematic contribution from the virial theorem. We aply this result 
for the COMA cluster data in next section. 



3 COMA Cluster Data Analysis 

The mass-to-light ratio of the COMA cluster was estimated, using the 
Newtonian theory, asTjvK 180 |5|. It seems reasonable to interpret this 
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large value for the mass-to-light ratio as a strong evidence in favour of 
dark matter. If MOND theory provide a succesfully description for the 
dynamics of this cluster we hope to find for T a value at least one 
order of magnitude less than 180. However, before any tentative to 
reconstruct a analysis in the COMA cluster using MOND theory, we 
need to assume some properties about the cluster: 

1 . We adopt a spherical symmetry. Although COMA have a elongated 
form, we are interested in the radial structure of the cluster. 

2. Any dark component of the system traces the luminous matter. 
With this assumption we consider a single component model for 
the cluster. This condition can be relaxed, however we will not 
consider this case. 

3. The cluster is in a stage of dynamical equilibrium and all the me- 
chanics of the cluster is well known and not yet relativistic. 

The above assumptions don't reflect actual gravitational systems, 
however will allow us to obtain in a simple way the cluster mass-to- 
light ratio. For a MOND description of orbits in non trivial symmetries 
see fT3l . 

The COMA cluster has been so widely studied over the years and 
for this reason there are a lot of available observational data in the 
literature. In order to obtain the velocity dispersion and the surface 
luminosity profile we take all the cluster members and binning them 
radially. We use the velocity dispersion data described in reference 
151 . Each bin contains aproximately 20 members of the cluster. The 
equivalent data for each bin and the best fit for the velocity dispersion 
and the surface luminosity are showed in figure 1. We do not adopt, 
a priori, any dynamical model for fitting the data. The best fit curves 
showed in figure 2 correspond to the curve that minimize the quatity \ 2 
that is proportional to the difference between the fitting value and the 
observational data in each bin. Inserting the surface luminosity profile 
obtained from the best fit, as we can see in figure 1, in the auxiliar 
expressions (16) we can compute the integral J, Eq. 1151 . With the 
velocity dispersion profile we can evaluate the kinetic term J in integral 
1121 . The mass-to-light ratio depends on the value of the critical radius 
r c . Table 1 shows results for Tm for various values of r c . 



4 Final Remarks 

The nature and origin of the dark matter is still one of the greatest chal- 
lenges in cosmology and astrophysics. Its influence in the dynamics of 
galaxies and clusters and its distribution in these systems has been 
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Figure 2: Best fit for the velocity dispersion and surface luminosity 
profiles. 

studied over the years and we don't have yet a satisfactory descrip- 
tion for this unknown component. In this contribution we apply the 
MOND theory in order to obtain a relation for the mass-to-light ratio 
in self- gravitating systems. Only with surface luminosity and disper- 
sion velocity data from COMA cluster the value of mass-to-light ratio 
was reduced using the MOND theory when compared with the Newto- 
nian one. The exact value of Tm depends on the critical radii r c . Table 
1 shows several results of our analysis. For suitable values of r c we 
found lowers mass-to-light ratios in contrast with the previous Newto- 
nian result sa 180|5|. The Newtonian limit of our relations (r c — > oo) 
indicates = 208.4 for the COMA cluster. We can admit that the 
estimation of the mass-to-light ratio as considered above is highly un- 
certain, however MOND correct the huge Newtonian mass-to-light ratio 
down to values compatible with baryonic matter alone. In particular, 



Table 1: Critical acceleration parameter (a ) and mass-to-light ratio 
in MOND theory (Tm)- We include the corresponding values for the 
critical radii r c , an auxiliar quantity derived from a . The Newtonian 
analysis (r c -> oo) reveals for COMA cluster a value for the mass-to- 
light ratio T N = 208.4 
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using the critical acceleration a ~ 10~ 10 m/,s 2 , a value that leads to a 
good fit of the rotation curve of spiral galaxies in the MOND theory |7|, 
a mass-to-light ratio of order of 70 is obtained, a third of the Newtonian 
value. In the analysis presented here, we have not taken into account 
the gas in the intergalactic medium of the cluster, which is typically re- 
sponsible for 20% of the total mass of a cluster 1141 . Even so, it agrees 
with the evaluation made in reference 1 15] where X-ray data are taken 
into account. Note, however, that our results disagree with the simpli- 
fied estimations made in |6| where just some general properties of the 
MOND theory have been used, leading to a very low mass/luminosity 
ratio, typically of order of a few tenths. Even if the necessity of dark 
matter remains, the dark matter problem is alleviate. Moreover, since 
there are many simplifications in the analysis done (the hypothesis of 
a continuous matter distribution, the choice of a spherical symmetry 
for the cluster, etc.), the results obtained using the interpolation func- 
tion (3) in the context of the MOND theory justify the investigation of a 
more realistic scenario, perhaps including a discrete distribution of the 
matter besides a continuous gas distribution. Of course, such more 
realistic analysis would imply the necessity of numerical simulations. 
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